X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;ds=sidebyside;f=PWGHF%2FvertexingHF%2FConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C;h=23558285f9162eac3cd617fae0d8036da7f3051d;hb=98667ce038b18d0e9f79bea854aed65558732c84;hp=1952efa3fa08eeec2c3534d782d1f8c1871d1c2f;hpb=2ac8110deab3cc81bd38a85ad0e2af52cba28aec;p=u%2Fmrichter%2FAliRoot.git diff --git a/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C b/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C index 1952efa3fa0..23558285f91 100644 --- a/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C +++ b/PWGHF/vertexingHF/ConfigVertexingHF_Pb_AllCent_NoLS_NoDstar_16var.C @@ -2,7 +2,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { printf("Call to AliAnalysisVertexingHF parameters setting :\n"); vHF = new AliAnalysisVertexingHF(); - + //--- switch-off candidates finding (default: all on) //vHF->SetD0toKpiOff(); vHF->SetJPSItoEleOff(); @@ -17,12 +17,13 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { vHF->SetCascadesOff(); vHF->SetFindVertexForCascades(kFALSE); vHF->SetMassCutBeforeVertexing(kTRUE); // PbPb + vHF->SetV0TypeForCascadeVertex(AliRDHFCuts::kAllV0s); - //--- set cuts for single-track selection - // displaced tracks + //--- set cuts for single-track selection + // displaced tracks AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default"); esdTrackCuts->SetRequireTPCRefit(kTRUE); - esdTrackCuts->SetMinNClustersTPC(70); + esdTrackCuts->SetMinNClustersTPC(50); esdTrackCuts->SetRequireITSRefit(kTRUE); //esdTrackCuts->SetMinNClustersITS(4); esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, @@ -40,14 +41,14 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { // displaced tracks for 20% mot central events 2 prongs AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default"); esdTrackCuts->SetRequireTPCRefit(kTRUE); - esdTrackCuts->SetMinNClustersTPC(70); + esdTrackCuts->SetMinNClustersTPC(50); esdTrackCuts->SetRequireITSRefit(kTRUE); //esdTrackCuts->SetMinNClustersITS(4); esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny); // |d0|>75 micron for pt<2GeV, no cut above 2 esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0050*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))"); - esdTrackCuts->SetMaxDCAToVertexXY(1.); + esdTrackCuts->SetMaxDCAToVertexXY(1.); esdTrackCuts->SetMaxDCAToVertexZ(1.); esdTrackCuts->SetPtRange(0.5,1.e10); esdTrackCuts->SetEtaRange(-0.8,+0.8); @@ -58,14 +59,14 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { // displaced tracks for 20% mot central events 3 prongs AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default"); esdTrackCuts->SetRequireTPCRefit(kTRUE); - esdTrackCuts->SetMinNClustersTPC(70); + esdTrackCuts->SetMinNClustersTPC(50); esdTrackCuts->SetRequireITSRefit(kTRUE); //esdTrackCuts->SetMinNClustersITS(4); esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny); // |d0|>75 micron for pt<2GeV, no cut above 2 esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0075*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))"); - esdTrackCuts->SetMaxDCAToVertexXY(1.); + esdTrackCuts->SetMaxDCAToVertexXY(1.); esdTrackCuts->SetMaxDCAToVertexZ(1.); esdTrackCuts->SetPtRange(0.6,1.e10); esdTrackCuts->SetEtaRange(-0.8,+0.8); @@ -76,7 +77,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { // D* soft pion tracks AliESDtrackCuts *esdTrackCutsSoftPi = new AliESDtrackCuts("AliESDtrackCuts","default"); esdTrackCutsSoftPi->SetMinNClustersITS(3); - esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.); + esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.); esdTrackCutsSoftPi->SetMaxDCAToVertexZ(1.); esdTrackCutsSoftPi->SetPtRange(0.1,1.e10); esdTrackCutsSoftPi->SetEtaRange(-0.8,+0.8); @@ -91,7 +92,8 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { cutsD0toKpi->SetUsePID(kFALSE); cutsD0toKpi->SetUseSpecialCuts(kFALSE); cutsD0toKpi->SetUsePhysicsSelection(kFALSE); - cutsD0toKpi->SetMaxCentrality(80.); + cutsD0toKpi->SetMinCentrality(-10.); + cutsD0toKpi->SetMaxCentrality(110.); cutsD0toKpi->SetMaxVtxZ(1.e6); cutsD0toKpi->SetTriggerClass(""); Float_t cutsArrayD0toKpi[11]={0.4,999999.,1.1,0.,0.,999999.,999999.,0.,0.5,-1,0.}; @@ -144,15 +146,15 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { cutsDStartoKpipi->SetUsePID(kFALSE); const Int_t nvars=16; const Int_t nptbins=2; - + Float_t* ptbins; ptbins=new Float_t[nptbins+1]; ptbins[0]=0.; ptbins[1]=5.; ptbins[2]=999.; - + cutsDStartoKpipi->SetPtBins(nptbins+1,ptbins); - + Float_t** rdcutsvalmine; rdcutsvalmine=new Float_t*[nvars]; for(Int_t iv=0;ivSetCuts(nvars,nptbins,rdcutsvalmine); - + cutsDStartoKpipi->AddTrackCuts(esdTrackCuts); cutsDStartoKpipi->AddTrackCutsSoftPi(esdTrackCutsSoftPi); cutsDStartoKpipi->SetMinPtCandidate(1.); @@ -203,16 +205,27 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { //-------------------------------------------------------- AliRDHFCutsLctoV0 *cutsLctoV0 = new AliRDHFCutsLctoV0("CutsLctoV0"); - Float_t cutsArrayLctoV0[9]={4.0,4.0,2.0,2.0,0.0,0.0,0.0,1000.,1000.}; - cutsLctoV0->SetCuts(9,cutsArrayLctoV0); + Float_t cutsArrayLctoV0[17]={1.0,1.0,0.05,0.05,0.0,0.0,0.0,1000.,1000.,0.99,3.,1000.,0.,0.,0.,0.,0.0}; + cutsLctoV0->SetCuts(17,cutsArrayLctoV0); cutsLctoV0->AddTrackCuts(esdTrackCuts); + + AliESDtrackCuts *esdV0daughterTrackCuts = new AliESDtrackCuts("AliESDtrackCutsForV0D","default cuts for V0 daughters"); + esdV0daughterTrackCuts->SetRequireTPCRefit(kTRUE); + esdV0daughterTrackCuts->SetMinNClustersTPC(30); + esdV0daughterTrackCuts->SetRequireITSRefit(kFALSE); + esdV0daughterTrackCuts->SetMinDCAToVertexXY(0.); + esdV0daughterTrackCuts->SetPtRange(0.05,1.e10); + esdV0daughterTrackCuts->SetEtaRange(-1.1,+1.1); + esdV0daughterTrackCuts->SetAcceptKinkDaughters(kTRUE); + esdV0daughterTrackCuts->SetRequireSigmaToVertex(kFALSE); + cutsLctoV0->AddTrackCutsV0daughters(esdV0daughterTrackCuts); vHF->SetCutsLctoV0(cutsLctoV0); - // + // //--- set this if you want to reconstruct primary vertex candidate by - // candidate using other tracks in the event (for pp, broad + // candidate using other tracks in the event (for pp, broad // interaction region) //vHF->SetRecoPrimVtxSkippingTrks(); - //--- OR set this if you want to remove the candidate daughters from + //--- OR set this if you want to remove the candidate daughters from // the primary vertex, without recostructing it from scratch //vHF->SetRmTrksFromPrimVtx(); @@ -221,7 +234,7 @@ AliAnalysisVertexingHF* ConfigVertexingHF() { //--- verbose // AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2); - + return vHF; }