]>
Commit | Line | Data |
---|---|---|
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->SetLikeSignOn(); // like-sign pairs and triplets | |
11 | vHF->Set4ProngOff(); | |
12 | //vHF->SetDstarOff(); | |
13 | vHF->SetFindVertexForDstar(kFALSE); | |
14 | //--- secondary vertex with KF? | |
15 | //vHF->SetSecVtxWithKF(); | |
16 | vHF->SetCascadesOff(); | |
17 | vHF->SetFindVertexForCascades(kFALSE); | |
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(4); | |
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 | Float_t cutsArrayD0toKpi[9]={0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.}; | |
44 | cutsD0toKpi->SetCuts(9,cutsArrayD0toKpi); | |
45 | cutsD0toKpi->AddTrackCuts(esdTrackCuts); | |
46 | vHF->SetCutsD0toKpi(cutsD0toKpi); | |
47 | AliRDHFCutsJpsitoee *cutsJpsitoee = new AliRDHFCutsJpsitoee("CutsJpsitoee"); | |
48 | Float_t cutsArrayJpsitoee[9]={0.350,100000.,1.1,0.,0.,100000.,100000.,100000000.,-1.1}; | |
49 | cutsJpsitoee->SetCuts(9,cutsArrayJpsitoee); | |
50 | cutsJpsitoee->AddTrackCuts(esdTrackCuts); | |
51 | vHF->SetCutsJpsitoee(cutsJpsitoee); | |
52 | AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi("CutsDplustoKpipi"); | |
53 | Float_t cutsArrayDplustoKpipi[12]={0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85,0.,10000000000.}; | |
54 | cutsDplustoKpipi->SetCuts(12,cutsArrayDplustoKpipi); | |
55 | cutsDplustoKpipi->AddTrackCuts(esdTrackCuts); | |
56 | vHF->SetCutsDplustoKpipi(cutsDplustoKpipi); | |
57 | AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi("CutsDstoKKpi"); | |
58 | Float_t cutsArrayDstoKKpi[14]={0.2,0.4,0.4,0.,0.,0.005,0.06,0.,0.,0.85,0.,1000,0.1,0.1}; | |
59 | cutsDstoKKpi->SetCuts(14,cutsArrayDstoKKpi); | |
60 | cutsDstoKKpi->AddTrackCuts(esdTrackCuts); | |
61 | vHF->SetCutsDstoKKpi(cutsDstoKKpi); | |
62 | AliRDHFCutsLctopKpi *cutsLctopKpi = new AliRDHFCutsLctopKpi("CutsLctopKpi"); | |
63 | Float_t cutsArrayLctopKpi[12]={0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85,0.,10000000000.}; | |
64 | cutsLctopKpi->SetCuts(12,cutsArrayLctopKpi); | |
65 | cutsLctopKpi->AddTrackCuts(esdTrackCuts); | |
66 | vHF->SetCutsLctopKpi(cutsLctopKpi); | |
67 | AliRDHFCutsD0toKpipipi *cutsD0toKpipipi = new AliRDHFCutsD0toKpipipi("CutsD0toKpipipi"); | |
68 | Float_t cutsArrayD0toKpipipi[9]={0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.}; | |
69 | cutsD0toKpipipi->SetCuts(9,cutsArrayD0toKpipipi); | |
70 | cutsD0toKpipipi->AddTrackCuts(esdTrackCuts); | |
71 | vHF->SetCutsD0toKpipipi(cutsD0toKpipipi); | |
72 | AliRDHFCutsD0toKpi *cutsD0fromDstar = new AliRDHFCutsD0toKpi("CutsD0fromDstar"); | |
73 | Float_t cutsArrayD0fromDstar[9]={0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.}; | |
74 | cutsD0fromDstar->SetCuts(9,cutsArrayD0fromDstar); | |
75 | cutsD0fromDstar->AddTrackCuts(esdTrackCuts); | |
76 | vHF->SetCutsD0fromDstar(cutsD0fromDstar); | |
77 | AliRDHFCutsLctoV0 *cutsLctoV0 = new AliRDHFCutsLctoV0("CutsLctoV0"); | |
78 | Float_t cutsArrayLctoV0[9]={4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000.,1000.}; | |
79 | cutsLctoV0->SetCuts(9,cutsArrayLctoV0); | |
80 | cutsLctoV0->AddTrackCuts(esdTrackCuts); | |
81 | vHF->SetCutsLctoV0(cutsLctoV0); | |
82 | // | |
83 | // to be removed: | |
84 | vHF->SetD0toKpiCuts(0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.); | |
85 | vHF->SetBtoJPSICuts(0.350,999999.,1.1,0.,0.,999999.,999999.,999999.,0.); | |
86 | vHF->SetDplusCuts(0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85,0,1e6); | |
87 | vHF->SetDsCuts(0.2,0.4,0.4,0.,0.,0.005,0.06,0.,0.,0.85,0.,0.1,0.1,0.1); | |
88 | vHF->SetLcCuts(0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85,0,1e6); | |
89 | vHF->SetD0to4ProngsCuts(0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.); | |
90 | vHF->SetDstarCuts(0.3, 0.1, 0.05, 100000000000.0, 0.5); | |
91 | vHF->SetD0fromDstarCuts(0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.); | |
92 | vHF->SetLctoV0Cuts(4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000,1000); | |
93 | //--- set this if you want to reconstruct primary vertex candidate by | |
94 | // candidate using other tracks in the event (for pp, broad | |
95 | // interaction region) | |
96 | //vHF->SetRecoPrimVtxSkippingTrks(); | |
97 | //--- OR set this if you want to remove the candidate daughters from | |
98 | // the primary vertex, without recostructing it from scratch | |
99 | //vHF->SetRmTrksFromPrimVtx(); | |
100 | ||
101 | //--- check the settings | |
102 | vHF->PrintStatus(); | |
103 | //--- verbose | |
104 | // AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2); | |
105 | ||
106 | ||
107 | return vHF; | |
108 | } | |
109 | ||
110 |