Int_t abspdgGranma =0;
Bool_t isFromB=kFALSE;
Bool_t isQuarkFound=kFALSE;
- while (mother >0 ){
+ while (mother >=0 ){
istep++;
AliDebug(2,Form("mother at step %d = %d", istep, mother));
AliAODMCParticle* mcGranma = dynamic_cast<AliAODMCParticle*>(fmcArray->At(mother));
Int_t label0 = fmcPartCandidate->GetDaughter(0);
Int_t label1 = fmcPartCandidate->GetDaughter(1);
AliDebug(3,Form("label0 = %d, label1 = %d",label0,label1));
- if (label1==0 || label0 == 0){
+ if (label1<=0 || label0 <= 0){
AliDebug(2, Form("The MC particle doesn't have correct daughters, skipping!!"));
return checkDaughters;
}
AliAODMCParticle *mcPartDaughter;
Int_t label0 = fmcPartCandidate->GetDaughter(0);
Int_t label1 = fmcPartCandidate->GetDaughter(1);
- if (label1==0 || label0 == 0){
+ if (label1<=0 || label0 <= 0){
AliDebug(2, Form("The MC particle doesn't have correct daughters, skipping!!"));
return bMCAccStep;
}
Int_t label0 = fmcPartCandidate->GetDaughter(0);
Int_t label1 = fmcPartCandidate->GetDaughter(1);
- if (label1==0 || label0 == 0){
+ if (label1<=0 || label0 <= 0){
AliDebug(2, Form("The MC particle doesn't have correct daughters, skipping!!"));
return bRefitStep;
}
}
//___________________________________________________________
-Bool_t AliCFVertexingHF::FillUnfoldingMatrix(Double_t fill[4]) const
+Bool_t AliCFVertexingHF::FillUnfoldingMatrix(UInt_t pdg, Double_t fill[4]) const
{
//
// filling the unfolding matrix
if(fmcPartCandidate){
fill[0] = GetPtCand();
- fill[1] = GetYCand();
+ fill[1] = GetYCand(pdg);
fill[2] = fmcPartCandidate->Pt();
fill[3] = fmcPartCandidate->Y();
Int_t label0 = fmcPartCandidate->GetDaughter(0);
Int_t label1 = fmcPartCandidate->GetDaughter(1);
AliDebug(2,Form("label0 = %d, label1 = %d",label0,label1));
- if (label1==0 || label0 == 0){
+ if (label1<=0 || label0 <= 0){
AliDebug(2, Form("The MC particle doesn't have correct daughters, skipping!!"));
delete [] fLabelArray;
fLabelArray = 0x0;