+ }
+
+ else if (shunt == 2) {
+ // the "primary" particle associated to the hit is
+ // the last track that has been flagged in the StepManager
+ // used by PHOS to associate the hit with the decay gamma
+ // rather than with the original pi0
+ TParticle *part;
+ Int_t current;
+ Int_t parent=track;
+ while (1) {
+ current=parent;
+ part = gAlice->GetMCApp()->Particle(current);
+ parent=part->GetFirstMother();
+ if(parent<0 || part->TestBit(kKeepBit))
+ break;
+ }
+ fTrack=current;
+ }
+
+ else {