Added counter for fake daughters

diff --git a/PWG3/vertexingHF/AliAODRecoDecayHF.h b/PWG3/vertexingHF/AliAODRecoDecayHF.h
index 5afe5d9e3acdd76e28e21a18c34f7729cd2aeb11..7ba02e3134d653f84881b290e62a23c5ce7356b2 100644 (file)
--- a/PWG3/vertexingHF/AliAODRecoDecayHF.h
+++ b/PWG3/vertexingHF/AliAODRecoDecayHF.h
@@ -114,6 +114,8 @@ class AliAODRecoDecayHF : public AliAODRecoDecay {
   Bool_t  HasSelectionBit(Int_t i) const {return TESTBIT(fSelectionMap,i);}
   ULong_t GetSelectionMap() const {return fSelectionMap;}
 
+  Int_t   NumberOfFakeDaughters() const;
+
  protected:
 
   AliAODVertex *fOwnPrimaryVtx; // primary vertex for this candidate
@@ -182,5 +184,27 @@ inline Bool_t AliAODRecoDecayHF::DaughterHasPointOnITSLayer(Int_t dg,Int_t l) co
   return TESTBIT(t->GetITSClusterMap(),l);
 }
 
+inline Int_t AliAODRecoDecayHF::NumberOfFakeDaughters() const 
+{
+  // Count number of daughters with negative label
+
+  Int_t nfakes=0;
+  for(Int_t i=0; i<GetNDaughters(); i++) {
+    AliAODTrack *track=(AliAODTrack*)GetDaughter(i);
+
+    if(track->Charge()==0) { // this is a two prong decay
+      AliAODRecoDecay *rd=(AliAODRecoDecay*)GetDaughter(i);
+      for(Int_t j=0; j<rd->GetNDaughters(); j++) {
+       AliAODTrack *track2=(AliAODTrack*)GetDaughter(j);
+       if(track2->GetLabel()<0) nfakes++;
+      }
+      continue;
+    }
+
+    if(track->GetLabel()<0) nfakes++;
+  }
+  return nfakes;
+}
+
 #endif