1 AliAnalysisVertexingHF* ConfigVertexingHF() {
3 printf("MYCONFIGPID Call to AliAnalysisVertexingHF parameters setting :\n");
4 vHF = new AliAnalysisVertexingHF();
6 //--- switch-off candidates finding (default: all on)
7 //vHF->SetD0toKpiOff();
9 vHF->SetJPSItoEleOff();
10 //vHF->Set3ProngOff();
11 //vHF->SetLikeSignOn(); // like-sign pairs and triplets
12 // vHF->SetLikeSign3prongOff();
15 vHF->SetFindVertexForDstar(kFALSE);
16 //--- secondary vertex with KF?
17 //vHF->SetSecVtxWithKF();
20 vHF->SetFindVertexForCascades(kTRUE);
21 vHF->SetV0TypeForCascadeVertex(AliRDHFCuts::kOnlyOfflineV0s); //All V0s 0, Offline 1, OnTheFly 2
22 vHF->SetUseProtonPIDforLambdaC2V0();
24 vHF->SetMassCutBeforeVertexing(kTRUE); // PbPb
27 vHF->SetUseKaonPIDfor3Prong(kTRUE);
28 vHF->SetUseProtonAndPionPIDforLambdaC();
29 vHF->SetnSigmaTOFforKaonSel(3., 5.);
30 vHF->SetnSigmaTPCforKaonSel(5., 5.);
31 vHF->SetnSigmaTOFforProtonSel(3.,5.);
32 vHF->SetnSigmaTPCforProtonSel(5., 5.);
33 vHF->SetnSigmaTPCforPionSel(4., 4.);
34 vHF->SetnSigmaTOFforPionSel(40.,40.);
35 vHF->SetMaxMomForTPCPid(9999999999.);
36 vHF->SetUseTPCPID(kTRUE);
37 vHF->SetUseTOFPID(kTRUE);
38 vHF->SetUseTPCPIDOnlyIfNoTOF(kTRUE);
40 //--- set cuts for single-track selection
42 AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
43 esdTrackCuts->SetRequireTPCRefit(kTRUE);
44 esdTrackCuts->SetMinNClustersTPC(70);
45 esdTrackCuts->SetRequireITSRefit(kTRUE);
46 //esdTrackCuts->SetMinNClustersITS(4);
47 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
48 AliESDtrackCuts::kAny);
49 // |d0|>30 micron for pt<2GeV, no cut above 2
50 esdTrackCuts->SetMinDCAToVertexXYPtDep("0.003*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))");
51 esdTrackCuts->SetMaxDCAToVertexXY(1.);
52 esdTrackCuts->SetMaxDCAToVertexZ(1.);
53 esdTrackCuts->SetPtRange(0.4,1.e10);
54 esdTrackCuts->SetEtaRange(-0.8,+0.8);
55 AliAnalysisFilter *trkFilter = new AliAnalysisFilter("trackFilter");
56 trkFilter->AddCuts(esdTrackCuts);
57 vHF->SetTrackFilter(trkFilter);
58 // D* soft pion tracks
59 AliESDtrackCuts *esdTrackCutsSoftPi = new AliESDtrackCuts("AliESDtrackCuts","default");
60 esdTrackCutsSoftPi->SetMinNClustersITS(4);
61 esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.);
62 esdTrackCutsSoftPi->SetMaxDCAToVertexZ(1.);
63 esdTrackCutsSoftPi->SetPtRange(0.2,1.e10);
64 esdTrackCutsSoftPi->SetEtaRange(-0.8,+0.8);
65 AliAnalysisFilter *trkFilterSoftPi = new AliAnalysisFilter("trackFilterSoftPi");
66 trkFilterSoftPi->AddCuts(esdTrackCutsSoftPi);
67 vHF->SetTrackFilterSoftPi(trkFilterSoftPi);
68 //--- set cuts for candidates selection
69 Int_t nptbins=2; Float_t ptlimits[2]={0.,1000000.};
70 AliRDHFCutsD0toKpi *cutsD0toKpi = new AliRDHFCutsD0toKpi("CutsD0toKpi");
71 cutsD0toKpi->SetStandardCutsPbPb2010();
72 cutsD0toKpi->SetMinCentrality(-10);
73 cutsD0toKpi->SetMaxCentrality(110);
74 cutsD0toKpi->SetUseSpecialCuts(kFALSE);
75 cutsD0toKpi->SetMinPtCandidate(0.);
76 cutsD0toKpi->SetUsePID(kFALSE);
77 cutsD0toKpi->SetUsePhysicsSelection(kFALSE);
78 cutsD0toKpi->SetMaxVtxZ(1.e6);
79 cutsD0toKpi->SetTriggerClass("");
80 Float_t cutsArrayD0toKpi[11]={0.4,999999.,1.1,0.,0.,999999.,999999.,0.,0.5,-1,0.};
81 cutsD0toKpi->SetPtBins(nptbins,ptlimits);
82 cutsD0toKpi->SetCuts(11,cutsArrayD0toKpi);
83 cutsD0toKpi->AddTrackCuts(esdTrackCuts);
84 vHF->SetCutsD0toKpi(cutsD0toKpi);
85 AliRDHFCutsJpsitoee *cutsJpsitoee = new AliRDHFCutsJpsitoee("CutsJpsitoee");
86 Float_t cutsArrayJpsitoee[9]={0.350,100000.,1.1,0.,0.,100000.,100000.,100000000.,-1.1};
87 cutsJpsitoee->SetCuts(9,cutsArrayJpsitoee);
88 cutsJpsitoee->AddTrackCuts(esdTrackCuts);
89 vHF->SetCutsJpsitoee(cutsJpsitoee);
90 AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi("CutsDplustoKpipi");
91 cutsDplustoKpipi->SetStandardCutsPbPb2010();
92 cutsDplustoKpipi->SetUsePID(kFALSE);
93 Float_t cutsArrayDplustoKpipi[14]={0.,0.3,0.3,0.,0.,0.01,0.05,0.05,0.,0.88,0.,10000000000.,0.,-1.};
94 cutsDplustoKpipi->SetPtBins(nptbins,ptlimits);
95 cutsDplustoKpipi->SetCuts(14,cutsArrayDplustoKpipi);
96 cutsDplustoKpipi->AddTrackCuts(esdTrackCuts);
97 cutsDplustoKpipi->SetMinPtCandidate(2000000000.);
98 vHF->SetCutsDplustoKpipi(cutsDplustoKpipi);
99 AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi("CutsDstoKKpi");
100 cutsDstoKKpi->SetStandardCutsPbPb2010();
101 cutsDstoKKpi->SetUsePID(kFALSE);
102 Float_t cutsArrayDstoKKpi[20]={0.,0.3,0.3,0.,0.,0.005,0.06,0.,0.,0.9,0.,100000.,0.035,0.0001,-1.,1.,0.,0.,0.,-1.};
103 cutsDstoKKpi->SetPtBins(nptbins,ptlimits);
104 cutsDstoKKpi->SetCuts(20,cutsArrayDstoKKpi);
105 cutsDstoKKpi->AddTrackCuts(esdTrackCuts);
106 cutsDstoKKpi->SetMinPtCandidate(1000000000.);
107 vHF->SetCutsDstoKKpi(cutsDstoKKpi);
108 AliRDHFCutsLctopKpi *cutsLctopKpi = new AliRDHFCutsLctopKpi("CutsLctopKpi");
109 cutsLctopKpi->SetStandardCutsPbPb2010();
110 cutsLctopKpi->SetUsePID(kFALSE);
111 Float_t cutsArrayLctopKpi[13]={0.13,0.4,0.4,0.,0.,0.,0.06,0.,0.,0.,0.,0.05,0.4};
112 cutsLctopKpi->SetPtBins(nptbins,ptlimits);
113 cutsLctopKpi->SetCuts(13,cutsArrayLctopKpi);
114 cutsLctopKpi->AddTrackCuts(esdTrackCuts);
115 cutsLctopKpi->SetMinPtCandidate(2.);
116 vHF->SetCutsLctopKpi(cutsLctopKpi);
117 AliRDHFCutsD0toKpipipi *cutsD0toKpipipi = new AliRDHFCutsD0toKpipipi("CutsD0toKpipipi");
118 Float_t cutsArrayD0toKpipipi[9]={0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.};
119 cutsD0toKpipipi->SetCuts(9,cutsArrayD0toKpipipi);
120 cutsD0toKpipipi->AddTrackCuts(esdTrackCuts);
121 vHF->SetCutsD0toKpipipi(cutsD0toKpipipi);
124 // D* pt dependent cuts ------------------------------------------
126 AliRDHFCutsDStartoKpipi *cutsDStartoKpipi = new AliRDHFCutsDStartoKpipi("CutsDStartoKpipi");
127 cutsDStartoKpipi->SetUsePID(kFALSE);
129 const Int_t nvars=16;
130 const Int_t nptbins=2;
133 ptbins=new Float_t[nptbins+1];
138 cutsDStartoKpipi->SetPtBins(nptbins+1,ptbins);
140 Float_t** rdcutsvalmine;
141 rdcutsvalmine=new Float_t*[nvars];
142 for(Int_t iv=0;iv<nvars;iv++){
143 rdcutsvalmine[iv]=new Float_t[nptbins];
146 rdcutsvalmine[0][0]=0.10; //D0 inv mass window
147 rdcutsvalmine[1][0]=0.06; // dca
148 rdcutsvalmine[2][0]=0.9; // thetastar
149 rdcutsvalmine[3][0]=0.5; // pt Pion
150 rdcutsvalmine[4][0]=0.5; // Pt Kaon
151 rdcutsvalmine[5][0]=0.1; // d0K
152 rdcutsvalmine[6][0]=0.1; // d0Pi
153 rdcutsvalmine[7][0]=0.0001; // d0xd0
154 rdcutsvalmine[8][0]=0.8; // costhetapoint
155 rdcutsvalmine[9][0]=0.15; // Dstar inv mass window
156 rdcutsvalmine[10][0]=0.03; // half width of (M_Kpipi-M_D0)
157 rdcutsvalmine[11][0]=0.1; // Pt min of Pi soft
158 rdcutsvalmine[12][0]=100.; // Pt max of pi soft
159 rdcutsvalmine[13][0]=9999.; // theta
160 rdcutsvalmine[14][0]=0.9; // |cosThetaPointXY|
161 rdcutsvalmine[15][0]=1.; // NormDecayLenghtXY
163 rdcutsvalmine[0][1]=0.10; //D0 inv mass window
164 rdcutsvalmine[1][1]=0.06; // dca
165 rdcutsvalmine[2][1]=0.9; // thetastar
166 rdcutsvalmine[3][1]=0.5; // pt Pion
167 rdcutsvalmine[4][1]=0.5; // Pt Kaon
168 rdcutsvalmine[5][1]=0.1; // d0K
169 rdcutsvalmine[6][1]=0.1; // d0Pi
170 rdcutsvalmine[7][1]=0.0001; // d0xd0
171 rdcutsvalmine[8][1]=0.7; // costhetapoint
172 rdcutsvalmine[9][1]=0.15; // Dstar inv mass window
173 rdcutsvalmine[10][1]=0.03; // half width of (M_Kpipi-M_D0)
174 rdcutsvalmine[11][1]=0.1; // Pt min of Pi soft
175 rdcutsvalmine[12][1]=100.; // Pt max of pi soft
176 rdcutsvalmine[13][1]=9999.; // theta
177 rdcutsvalmine[14][1]=0.8; // |cosThetaPointXY|
178 rdcutsvalmine[15][1]=0.; // NormDecayLenghtXY
180 cutsDStartoKpipi->SetCuts(nvars,nptbins,rdcutsvalmine);
182 cutsDStartoKpipi->AddTrackCuts(esdTrackCuts);
183 cutsDStartoKpipi->AddTrackCutsSoftPi(esdTrackCutsSoftPi);
184 cutsDStartoKpipi->SetMinPtCandidate(2.);
185 vHF->SetCutsDStartoKpipi(cutsDStartoKpipi);
187 //--------------------------------------------------------
189 AliRDHFCutsLctoV0 *cutsLctoV0 = new AliRDHFCutsLctoV0("CutsLctoV0");
191 Float_t cutsArrayLctoV0[9]={4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000.,1000.};
192 cutsLctoV0->SetCuts(9,cutsArrayLctoV0);
195 Float_t cutsArrayLctoV0[10]={4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000.,1000.,1.};
196 cutsLctoV0->SetCuts(10,cutsArrayLctoV0);
197 cutsLctoV0->AddTrackCuts(esdTrackCuts);
198 vHF->SetCutsLctoV0(cutsLctoV0);
200 //--- set this if you want to reconstruct primary vertex candidate by
201 // candidate using other tracks in the event (for pp, broad
202 // interaction region)
203 //vHF->SetRecoPrimVtxSkippingTrks();
204 //--- OR set this if you want to remove the candidate daughters from
205 // the primary vertex, without recostructing it from scratch
206 //vHF->SetRmTrksFromPrimVtx();
208 //--- check the settings
211 // AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2);