Added cut |zVtxMC|<fRDHFCuts->GetMaxVtxZ() (Chiara)

diff --git a/PWG3/vertexingHF/AliCFTaskVertexingHF.cxx b/PWG3/vertexingHF/AliCFTaskVertexingHF.cxx
index 22c0774c1ced5b17a0004deb0c55ec7318a4c21f..2074ea827b1d8e5b0fc1d852bdb001ceeb3a46da 100644 (file)
--- a/PWG3/vertexingHF/AliCFTaskVertexingHF.cxx
+++ b/PWG3/vertexingHF/AliCFTaskVertexingHF.cxx
@@ -451,6 +451,11 @@ void AliCFTaskVertexingHF::UserExec(Option_t *)
        Double_t zPrimVertex = aodVtx ->GetZ();
        Double_t zMCVertex = mcHeader->GetVtxZ();
        
+       if (TMath::Abs(zMCVertex) > fCuts->GetMaxVtxZ()){
+         AliDebug(3,Form("z coordinate of MC vertex = %f, it was required to be within [-%f, +%f], skipping event", zMCVertex, fCuts->GetMaxVtxZ(), fCuts->GetMaxVtxZ()));
+         return;
+       }
+
        AliESDtrackCuts** trackCuts = new AliESDtrackCuts*[cfVtxHF->GetNProngs()];
        if (fDecayChannel == 21){
                // for the D*, setting the third element of the array of the track cuts to those for the soft pion