1 AliAnalysisVertexingHF* ConfigVertexingHF() {
3 printf("Call to AliAnalysisVertexingHF parameters setting :\n");
4 vHF = new AliAnalysisVertexingHF();
6 //--- switch-off candidates finding (default: all on)
7 //vHF->SetD0toKpiOff();
8 vHF->SetJPSItoEleOff();
10 vHF->SetLikeSignOn(); // like-sign pairs and triplets
13 vHF->SetFindVertexForDstar(kFALSE);
14 //--- secondary vertex with KF?
15 //vHF->SetSecVtxWithKF();
16 vHF->SetCascadesOff();
17 vHF->SetFindVertexForCascades(kFALSE);
18 vHF->SetMassCutBeforeVertexing(kTRUE); // PbPb
20 //--- set cuts for single-track selection
22 AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
23 esdTrackCuts->SetRequireTPCRefit(kTRUE);
24 esdTrackCuts->SetMinNClustersTPC(70);
25 esdTrackCuts->SetRequireITSRefit(kTRUE);
26 //esdTrackCuts->SetMinNClustersITS(4);
27 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
28 AliESDtrackCuts::kAny);
29 // |d0|>100 micron for pt<2GeV, no cut above 2
30 esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0100*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))");
31 esdTrackCuts->SetMaxDCAToVertexXY(1.);
32 esdTrackCuts->SetMaxDCAToVertexZ(1.);
33 esdTrackCuts->SetPtRange(0.8,1.e10);
34 esdTrackCuts->SetEtaRange(-0.8,+0.8);
35 AliAnalysisFilter *trkFilter = new AliAnalysisFilter("trackFilter");
36 trkFilter->AddCuts(esdTrackCuts);
37 vHF->SetTrackFilter(trkFilter);
38 // D* soft pion tracks
39 AliESDtrackCuts *esdTrackCutsSoftPi = new AliESDtrackCuts("AliESDtrackCuts","default");
40 esdTrackCutsSoftPi->SetMinNClustersITS(4);
41 esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.);
42 esdTrackCutsSoftPi->SetMaxDCAToVertexZ(1.);
43 esdTrackCutsSoftPi->SetPtRange(0.2,1.e10);
44 esdTrackCutsSoftPi->SetEtaRange(-0.8,+0.8);
45 AliAnalysisFilter *trkFilterSoftPi = new AliAnalysisFilter("trackFilterSoftPi");
46 trkFilterSoftPi->AddCuts(esdTrackCutsSoftPi);
47 vHF->SetTrackFilterSoftPi(trkFilterSoftPi);
48 //--- set cuts for candidates selection
49 AliRDHFCutsD0toKpi *cutsD0toKpi = new AliRDHFCutsD0toKpi("CutsD0toKpi");
50 cutsD0toKpi->SetStandardCutsPbPb2010();
51 cutsD0toKpi->SetUsePID(kFALSE);
52 Float_t cutsArrayD0toKpi[9]={0.2,999999.,1.1,0.,0.,999999.,999999.,0.,0.5};
53 cutsD0toKpi->SetCuts(9,cutsArrayD0toKpi);
54 cutsD0toKpi->AddTrackCuts(esdTrackCuts);
55 vHF->SetCutsD0toKpi(cutsD0toKpi);
56 AliRDHFCutsJpsitoee *cutsJpsitoee = new AliRDHFCutsJpsitoee("CutsJpsitoee");
57 Float_t cutsArrayJpsitoee[9]={0.350,100000.,1.1,0.,0.,100000.,100000.,100000000.,-1.1};
58 cutsJpsitoee->SetCuts(9,cutsArrayJpsitoee);
59 cutsJpsitoee->AddTrackCuts(esdTrackCuts);
60 vHF->SetCutsJpsitoee(cutsJpsitoee);
61 AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi("CutsDplustoKpipi");
62 cutsDplustoKpipi->SetStandardCutsPbPb2010();
63 cutsDplustoKpipi->SetUsePID(kFALSE);
64 Float_t cutsArrayDplustoKpipi[12]={0.2,0.8,0.8,0.,0.,0.01,0.06,0.03,0.,0.85,0.,10000000000.};
65 cutsDplustoKpipi->SetCuts(12,cutsArrayDplustoKpipi);
66 cutsDplustoKpipi->AddTrackCuts(esdTrackCuts);
67 vHF->SetCutsDplustoKpipi(cutsDplustoKpipi);
68 AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi("CutsDstoKKpi");
69 cutsDstoKKpi->SetStandardCutsPbPb2010();
70 cutsDstoKKpi->SetUsePID(kFALSE);
71 Float_t cutsArrayDstoKKpi[14]={0.,999.,999.,0.,0.,0.005,0.06,0.,0.,0.7,0.,1000.,0.1,0.1};
72 cutsDstoKKpi->SetCuts(14,cutsArrayDstoKKpi);
73 cutsDstoKKpi->AddTrackCuts(esdTrackCuts);
74 vHF->SetCutsDstoKKpi(cutsDstoKKpi);
75 AliRDHFCutsLctopKpi *cutsLctopKpi = new AliRDHFCutsLctopKpi("CutsLctopKpi");
76 cutsLctopKpi->SetStandardCutsPbPb2010();
77 cutsLctopKpi->SetUsePID(kFALSE);
78 Float_t cutsArrayLctopKpi[12]={0.13,0.5,0.6,0.,0.,0.01,0.04,0.006,0.8,0.3,0.,0.05};
79 cutsLctopKpi->SetCuts(12,cutsArrayLctopKpi);
80 cutsLctopKpi->AddTrackCuts(esdTrackCuts);
81 vHF->SetCutsLctopKpi(cutsLctopKpi);
82 AliRDHFCutsD0toKpipipi *cutsD0toKpipipi = new AliRDHFCutsD0toKpipipi("CutsD0toKpipipi");
83 Float_t cutsArrayD0toKpipipi[9]={0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.};
84 cutsD0toKpipipi->SetCuts(9,cutsArrayD0toKpipipi);
85 cutsD0toKpipipi->AddTrackCuts(esdTrackCuts);
86 vHF->SetCutsD0toKpipipi(cutsD0toKpipipi);
87 AliRDHFCutsDStartoKpipi *cutsDStartoKpipi = new AliRDHFCutsDStartoKpipi("CutsDStartoKpipi");
88 cutsDStartoKpipi->SetStandardCutsPbPb2010();
89 cutsDStartoKpipi->SetUsePID(kFALSE);
90 Float_t cutsArrayDStartoKpipi[14]={0.15,0.07,0.85,0.8,0.8,0.06.,0.06.,0.001,0.6, 0.15, 0.03, 0.2, 5, 0.5}; // first 9 for D0 from D*, last 5 for D*
91 cutsDStartoKpipi->SetCuts(14,cutsArrayDStartoKpipi);
92 cutsDStartoKpipi->AddTrackCuts(esdTrackCuts);
93 cutsDStartoKpipi->AddTrackCutsSoftPi(esdTrackCutsSoftPi);
94 vHF->SetCutsDStartoKpipi(cutsDStartoKpipi);
95 AliRDHFCutsLctoV0 *cutsLctoV0 = new AliRDHFCutsLctoV0("CutsLctoV0");
96 Float_t cutsArrayLctoV0[9]={4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000.,1000.};
97 cutsLctoV0->SetCuts(9,cutsArrayLctoV0);
98 cutsLctoV0->AddTrackCuts(esdTrackCuts);
99 vHF->SetCutsLctoV0(cutsLctoV0);
101 //--- set this if you want to reconstruct primary vertex candidate by
102 // candidate using other tracks in the event (for pp, broad
103 // interaction region)
104 //vHF->SetRecoPrimVtxSkippingTrks();
105 //--- OR set this if you want to remove the candidate daughters from
106 // the primary vertex, without recostructing it from scratch
107 //vHF->SetRmTrksFromPrimVtx();
109 //--- check the settings
112 // AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2);