+ if( ((depositedEnergy = gMC->Edep()) > 0.) && (gMC->TrackTime() < fTimeCut)){// Track is inside a scintillator and deposits some energy
+ if (fCurPrimary==-1)
+ fCurPrimary=gAlice->GetMCApp()->GetPrimary(tracknumber);
+
+ if (fCurParent==-1 || tracknumber != fCurTrack) {
+ // Check parentage
+ //Int_t parent=tracknumber;
+ parent=tracknumber;
+ if (fCurParent != -1) {
+ while (parent != fCurParent && parent != -1) {
+ //TParticle *part=gAlice->GetMCApp()->Particle(parent);
+ part=gAlice->GetMCApp()->Particle(parent);
+ parent=part->GetFirstMother();
+ }
+ }
+ if (fCurParent==-1 || parent==-1) {
+ //Int_t parent=tracknumber;
+ //TParticle *part=gAlice->GetMCApp()->Particle(parent);
+ parent=tracknumber;
+ part=gAlice->GetMCApp()->Particle(parent);
+ while (parent != -1 && geom->IsInEMCAL(part->Vx(),part->Vy(),part->Vz())) {
+ parent=part->GetFirstMother();
+ if (parent!=-1)
+ part=gAlice->GetMCApp()->Particle(parent);
+ }
+ fCurParent=parent;
+ if (fCurParent==-1)
+ Error("StepManager","Cannot find parent");
+ else {
+ //TParticle *part=gAlice->GetMCApp()->Particle(fCurParent);
+ part=gAlice->GetMCApp()->Particle(fCurParent);
+ ienergy = part->Energy();
+ }
+ while (parent != -1) {
+ part=gAlice->GetMCApp()->Particle(parent);
+ part->SetBit(kKeepBit);
+ parent=part->GetFirstMother();
+ }
+ }
+ fCurTrack=tracknumber;
+ }