AliPhysics  5eaf189 (5eaf189)
ConfigVertexingHF_pPbRHF_DplustoK0s_DstoK0s.C
Go to the documentation of this file.
2  printf("Call to AliAnalysisVertexingHF parameters setting :\n");
3  vHF = new AliAnalysisVertexingHF();
4  vHF->SetMakeReducedRHF(kTRUE);
5  //--- switch-off candidates finding (default: all on)
6  //vHF->SetD0toKpiOff();
7  vHF->SetJPSItoEleOff();
8  //vHF->Set3ProngOff();
9 // vHF->SetLikeSignOff(); // like-sign pairs and triplets
10  vHF->Set4ProngOff();
11  //vHF->SetDstarOff();
12  vHF->SetFindVertexForDstar(kFALSE);
13  //--- secondary vertex with KF?
14  //vHF->SetSecVtxWithKF();
15 // vHF->SetCascadesOff();
16  vHF->SetFindVertexForCascades(kFALSE);
17  vHF->SetV0TypeForCascadeVertex(AliRDHFCuts::kAllV0s);
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(50);
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");
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);
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.4,0.5,0.,0.,0.01,0.06,0.005,0.7,0.,0.,0.05,0.4};
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");
83  Float_t cutsArrayLctoV0[21]={1.0,1.0,0.05,0.05,0.0,0.0,0.0,1000.,1000.,0.99,3.,1000.,0.,0.,0.,0.,9999.,-9999.,-9999.,-9999.,0.0};
84  cutsLctoV0->SetCuts(21,cutsArrayLctoV0);
85  cutsLctoV0->AddTrackCuts(esdTrackCuts);
86  AliESDtrackCuts *esdV0daughterTrackCuts = new AliESDtrackCuts("AliESDtrackCutsForV0D","default cuts for V0 daughters");
87  esdV0daughterTrackCuts->SetRequireTPCRefit(kTRUE);
88  esdV0daughterTrackCuts->SetMinNClustersTPC(30);
89  esdV0daughterTrackCuts->SetRequireITSRefit(kFALSE);
90  esdV0daughterTrackCuts->SetMinDCAToVertexXY(0.);
91  esdV0daughterTrackCuts->SetPtRange(0.05,1.e10);
92  esdV0daughterTrackCuts->SetEtaRange(-1.1,+1.1);
93  esdV0daughterTrackCuts->SetAcceptKinkDaughters(kTRUE);
94  esdV0daughterTrackCuts->SetRequireSigmaToVertex(kFALSE);
95  cutsLctoV0->AddTrackCutsV0daughters(esdV0daughterTrackCuts);
96  vHF->SetCutsLctoV0(cutsLctoV0);
97 
98  AliRDHFCutsDplustoK0spi *cutsDplustoK0s = new AliRDHFCutsDplustoK0spi("CutsDplustoK0spi");
99  Float_t cutsArrayDplustoK0s[16] = {0.3, 0.2, 0. , 5., 0., 0.99, 0., 0., 0., 3., 20., 1.1, 0., 0., 0., 0.};
100  cutsDplustoK0s->SetCuts(16, cutsArrayDplustoK0s);
101  cutsDplustoK0s->AddTrackCuts(esdTrackCuts);
102  cutsDplustoK0s->AddTrackCutsV0daughters(esdV0daughterTrackCuts);
103  vHF->SetCutsDplustoK0spi(cutsDplustoK0s);
104 
105  AliRDHFCutsDstoK0sK *cutsDstoK0s = new AliRDHFCutsDstoK0sK("CutsDstoK0sK");
106  Float_t cutsArrayDstoK0s[16] = {0.3, 0.2, 0. , 5., 0., 0.99, 0., 0., 0., 3., 20., 1.1, 0., 0., 0., 0.};
107  cutsDstoK0s->SetCuts(16, cutsArrayDstoK0s);
108  cutsDstoK0s->AddTrackCuts(esdTrackCuts);
109  cutsDstoK0s->AddTrackCutsV0daughters(esdV0daughterTrackCuts);
110  vHF->SetCutsDstoK0sK(cutsDstoK0s);
111 
112  //
113  //--- set this if you want to reconstruct primary vertex candidate by
114  // candidate using other tracks in the event (for pp, broad
115  // interaction region)
116  //vHF->SetRecoPrimVtxSkippingTrks();
117  //--- OR set this if you want to remove the candidate daughters from
118  // the primary vertex, without recostructing it from scratch
119  //vHF->SetRmTrksFromPrimVtx();
120 
121  //--- check the settings
122  vHF->PrintStatus();
123  //--- verbose
124  // AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2);
125 
126 
127  return vHF;
128 }
129 
130 
void SetMaxVtxZ(Float_t z=1e6)
Definition: AliRDHFCuts.h:61
void AddTrackCutsV0daughters(AliESDtrackCuts *v0daug)
AliAnalysisVertexingHF * ConfigVertexingHF()
void SetUsePhysicsSelection(Bool_t use=kTRUE)
Definition: AliRDHFCuts.h:351
Class for cuts on AOD reconstructed D+->Kpipi.
void SetCuts(Int_t nVars, Int_t nPtBins, Float_t **cutsRD)
float Float_t
Definition: External.C:68
Class for cuts on AOD reconstructed Ds->K0S+K.
void AddTrackCutsV0daughters(AliESDtrackCuts *v0daug)
void AddTrackCutsSoftPi(const AliESDtrackCuts *cuts)
Class for cuts on AOD reconstructed D+->K0S+pi.
void AddTrackCutsV0daughters(AliESDtrackCuts *v0daug)
void AddTrackCuts(const AliESDtrackCuts *cuts)
Definition: AliRDHFCuts.h:207
void SetTriggerClass(TString trclass0, TString trclass1="")
Definition: AliRDHFCuts.h:197