Updates in PbPb cuts (Andrea Rossi)
[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliCFVertexingHF2Prong.cxx
index a0068bc..afdc999 100644 (file)
@@ -91,6 +91,16 @@ Bool_t AliCFVertexingHF2Prong::SetRecoCandidateParam(AliAODRecoDecayHF *recoCand
  
        Int_t mcLabel = fRecoCandidate->MatchToMC(pdgCand,fmcArray,2,pdgDgD0toKpi);
        if (mcLabel == -1) return bSignAssoc;
+
+       if (fRecoCandidate->NumberOfFakeDaughters()>0){
+               fFake = 0;    // fake candidate
+               if (fFakeSelection==1) return bSignAssoc;
+       }
+       if (fRecoCandidate->NumberOfFakeDaughters()==0){
+               fFake = 2;    // non-fake candidate
+               if (fFakeSelection==2) return bSignAssoc;
+       }
+
        SetMCLabel(mcLabel);
        fmcPartCandidate = dynamic_cast<AliAODMCParticle*>(fmcArray->At(fmcLabel));
        if (!fmcPartCandidate){
@@ -198,6 +208,7 @@ Bool_t AliCFVertexingHF2Prong::GetGeneratedValuesFromMCParticle(Double_t* vector
        vectorMC[11] = fmcPartCandidate->Phi(); 
        vectorMC[12] = fzMCVertex;    // z of reconstructed of primary vertex
        vectorMC[13] = fCentValue;   //reconstructed centrality 
+       vectorMC[14] = 1.;           // always filling with 1 at MC level 
 
        delete decay;
        bGenValues = kTRUE;
@@ -264,6 +275,7 @@ Bool_t AliCFVertexingHF2Prong::GetRecoValuesFromCandidate(Double_t *vectorReco)
        vectorReco[11] = phi;  
        vectorReco[12] = fzPrimVertex;    // z of reconstructed of primary vertex
        vectorReco[13] = fCentValue; //reconstructed centrality 
+       vectorReco[14] = fFake;      // whether the reconstructed candidate was a fake (fFake = 0) or not (fFake = 2) 
 
        bFillRecoValues = kTRUE;