From aed46c569aecbdd2f0eefad1d2b6b815a24deb9a Mon Sep 17 00:00:00 2001 From: dainese Date: Fri, 31 Jul 2009 11:11:51 +0000 Subject: [PATCH] Check kTPCrefit for daughter tracks and apply minimal D0 cuts common to LHC09a4 and LHC09a6 (as in ConfigVertexingHF.C) Chiara Z --- .../AliCFHeavyFlavourTaskMultiVarMultiStep.cxx | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) diff --git a/PWG3/vertexingHF/AliCFHeavyFlavourTaskMultiVarMultiStep.cxx b/PWG3/vertexingHF/AliCFHeavyFlavourTaskMultiVarMultiStep.cxx index b8a27ac296b..6b2788c2a81 100644 --- a/PWG3/vertexingHF/AliCFHeavyFlavourTaskMultiVarMultiStep.cxx +++ b/PWG3/vertexingHF/AliCFHeavyFlavourTaskMultiVarMultiStep.cxx @@ -46,6 +46,7 @@ #include "AliAODRecoDecayHF.h" #include "AliAODRecoDecayHF2Prong.h" #include "AliAODMCParticle.h" +#include "AliESDtrack.h" #include "TChain.h" #include "THnSparse.h" #include "TH2D.h" @@ -310,6 +311,21 @@ void AliCFHeavyFlavourTaskMultiVarMultiStep::UserExec(Option_t *) continue; } + // check whether the daughters have kTPCrefit set + AliAODTrack *track0 = (AliAODTrack*)d0tokpi->GetDaughter(0); + AliAODTrack *track1 = (AliAODTrack*)d0tokpi->GetDaughter(1); + if((!(track0->GetStatus()&AliESDtrack::kTPCrefit)) || (!(track1->GetStatus()&AliESDtrack::kTPCrefit))) { + // skipping if at least one daughter does not have kTPCrefit + continue; + } + + const Double_t d0tokpi_cuts[9] = {0.3,999999.,1.1,0.,0.,999999.,999999.,999999.,0.}; + Int_t okD0, okD0bar; + if (!(d0tokpi->SelectD0(&d0tokpi_cuts[0],okD0,okD0bar))){ + // skipping candidate + continue; + } + // check if associated MC v0 passes the cuts if (!fCFManager->CheckParticleCuts(0 ,mcVtxHF)) { // 0 stands for MC AliDebug(2, "Skipping the particles due to cuts"); -- 2.43.0