// build the combinatorics
Int_t nSelected=0;
Int_t nFiltered=0;
-
- AliAODVertex* v2=new AliAODVertex();
- AliAODVertex* v3=new AliAODVertex();
+ Double_t dummypos[3]={0.,0.,0.};
+ AliAODVertex* v2=new AliAODVertex(dummypos,999.,-1,2);
+ AliAODVertex* v3=new AliAODVertex(dummypos,999.,-1,3);
// dummy values of track impact parameter, needed to build an AliAODRecoDecay object
Double_t d02[2]={0.,0.};
Double_t d03[3]={0.,0.,0.};
Int_t pdgCode = TMath::Abs( part->GetPdgCode() );
if(pdgCode==321){
AliAODMCParticle* dmes = dynamic_cast<AliAODMCParticle*>(arrayMC->At(labD));
- if(dmes) fPtVsYReco->Fill(dmes->Pt(),dmes->Y());
+ if(dmes){
+ Int_t orig=AliVertexingHFUtils::CheckOrigin(arrayMC,dmes,fGoUpToQuark);
+ if((fPromptFeeddown==kFeeddown && orig==5)|| (fPromptFeeddown==kPrompt && orig==4) || (fPromptFeeddown==kBoth && orig>=4)) {
+ fPtVsYReco->Fill(dmes->Pt(),dmes->Y());
+ }
+ }
fMassVsPtVsYSig->Fill(mass,pt,rapid);
}else{
fMassVsPtVsYRefl->Fill(mass,pt,rapid);