- mcp->fIso = AliAnalysisTaskEMCALPhoton::GetMcIsolation(mcP, itrack, 0.4 , 0.2);
+ if(itrack == 8)
+ mcp->fIso = AliAnalysisTaskEMCALPhoton::GetMcIsolation(mcP, itrack, 0.4 , 0.2);
+ else
+ mcp->fIso = -1;
+ if(fDebug)
+ printf("\t::FillMcParts(): label=%d, pdg=%d, pt=%1.1f, mom=%d\n\t::FillMcParts() returning...\n\n", mcp->fLabel,mcp->fPdg,mcp->fPt,mcp->fMother);
+ for(Int_t id=mcP->GetFirstDaughter(); id <= mcP->GetLastDaughter(); id++){
+ if(id<=mcP->GetMother(0))
+ continue;
+ if(id<0 || id>nTracks)
+ continue;
+ TParticle *mcD = static_cast<TParticle*>(fStack->Particle(id));
+ if(!mcD)
+ continue;
+ FillMcPart(mcD, fMyMcIndex++, id);
+ }
+