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