New classes for the cuts on charm candidates
[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / ConfigVertexingHF_highmult.C
CommitLineData
289fa632 1AliAnalysisVertexingHF* ConfigVertexingHF() {
2f1e74b9 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
17 //--- set cuts for single-track selection
18 // displaced tracks
19 AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
20 esdTrackCuts->SetRequireTPCRefit(kTRUE);
21 esdTrackCuts->SetRequireITSRefit(kTRUE);
22 esdTrackCuts->SetMinNClustersITS(5);
23 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
24 AliESDtrackCuts::kBoth);
25 esdTrackCuts->SetMinDCAToVertexXY(0.);
26 esdTrackCuts->SetPtRange(0.3,1.e10);
27 AliAnalysisFilter *trkFilter = new AliAnalysisFilter("trackFilter");
28 trkFilter->AddCuts(esdTrackCuts);
29 vHF->SetTrackFilter(trkFilter);
30 // D* soft pion tracks
31 AliESDtrackCuts *esdTrackCutsSoftPi = new AliESDtrackCuts("AliESDtrackCuts","default");
32 esdTrackCutsSoftPi->SetMinNClustersITS(4);
33 AliAnalysisFilter *trkFilterSoftPi = new AliAnalysisFilter("trackFilterSoftPi");
34 trkFilterSoftPi->AddCuts(esdTrackCutsSoftPi);
35 vHF->SetTrackFilterSoftPi(trkFilterSoftPi);
36 //--- set cuts for candidates selection
e3d40058 37 AliRDHFCutsD0toKpi *cutsD0toKpi = new AliRDHFCutsD0toKpi();
38 Float_t cutsArrayD0toKpi[9]={0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.};
39 cutsD0toKpi->SetCuts(9,cutsArrayD0toKpi);
40 cutsD0toKpi->AddTrackCuts(esdTrackCuts);
41 vHF->SetCutsD0toKpi(cutsD0toKpi);
42 AliRDHFCutsJpsitoee *cutsJpsitoee = new AliRDHFCutsJpsitoee();
43 Float_t cutsArrayJpsitoee[9]={0.35,999999.,1.1,0.,0.,999999.,999999.,999999.,0.};
44 cutsJpsitoee->SetCuts(9,cutsArrayJpsitoee);
45 cutsJpsitoee->AddTrackCuts(esdTrackCuts);
46 vHF->SetCutsJpsitoee(cutsJpsitoee);
47 AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi();
48 Float_t cutsArrayDplustoKpipi[12]={0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85,0,1e6};
49 cutsDplustoKpipi->SetCuts(12,cutsArrayDplustoKpipi);
50 cutsDplustoKpipi->AddTrackCuts(esdTrackCuts);
51 vHF->SetCutsDplustoKpipi(cutsDplustoKpipi);
52 AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi();
53 Float_t cutsArrayDstoKKpi[14]={0.2,0.4,0.4,0.,0.,0.005,0.06,0.,0.,0.85,0.,0.1,0.1,0.1};
54 cutsDstoKKpi->SetCuts(14,cutsArrayDstoKKpi);
55 cutsDstoKKpi->AddTrackCuts(esdTrackCuts);
56 vHF->SetCutsDstoKKpi(cutsDstoKKpi);
57 AliRDHFCutsLctopKpi *cutsLctopKpi = new AliRDHFCutsLctopKpi();
58 Float_t cutsArrayLctopKpi[12]={0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85,0,1e6};
59 cutsLctopKpi->SetCuts(12,cutsArrayLctopKpi);
60 cutsLctopKpi->AddTrackCuts(esdTrackCuts);
61 vHF->SetCutsLctopKpi(cutsLctopKpi);
62 AliRDHFCutsD0toKpipipi *cutsD0toKpipipi = new AliRDHFCutsD0toKpipipi();
63 Float_t cutsArrayD0toKpipipi[9]={0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.};
64 cutsD0toKpipipi->SetCuts(9,cutsArrayD0toKpipipi);
65 cutsD0toKpipipi->AddTrackCuts(esdTrackCuts);
66 vHF->SetCutsD0toKpipipi(cutsD0toKpipipi);
67 AliRDHFCutsD0toKpi *cutsD0fromDstar = new AliRDHFCutsD0toKpi();
68 Float_t cutsArrayD0fromDstar[9]={0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.};
69 cutsD0fromDstar->SetCuts(9,cutsArrayD0fromDstar);
70 cutsD0fromDstar->AddTrackCuts(esdTrackCuts);
71 vHF->SetCutsD0fromDstar(cutsD0fromDstar);
72 //
2f1e74b9 73 vHF->SetD0toKpiCuts(0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.);
74 vHF->SetBtoJPSICuts(0.350);
75 vHF->SetDplusCuts(0.2,0.4,0.4,0.,0.,0.01,0.06,0.02,0.,0.85);
76 vHF->SetDsCuts(0.2,0.4,0.4,0.,0.,0.005,0.06,0.,0.,0.85,0.,0.1,0.1);
77 vHF->SetLcCuts(0.2,0.4,0.4,0.,0.,0.01,0.06,0.,0.,0.85);
78 vHF->SetD0to4ProngsCuts(0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.);
79 vHF->SetDstarCuts(0.3, 0.1, 0.05, 100000000000.0, 0.5);
80 vHF->SetD0fromDstarCuts(0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.);
81 //--- set this if you want to reconstruct primary vertex candidate by
82 // candidate using other tracks in the event (for pp, broad
83 // interaction region)
84 //vHF->SetRecoPrimVtxSkippingTrks();
85 //--- OR set this if you want to remove the candidate daughters from
86 // the primary vertex, without recostructing it from scratch
87 //vHF->SetRmTrksFromPrimVtx();
88
89 //--- check the settings
90 vHF->PrintStatus();
91 //--- verbose
92 //AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2);
93
94
95 return vHF;
96}
97
98
99